The Corley-Jacobson dispersion relation and trans- Planckian inflation 
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In this Letter we study the dependence of the spectrum of fluctuations in inflationary cosmology 
on possible effects of trans-Planckian physics, using the Corley/Jacobson dispersion relations as 
an example. We compare the methods used in previous work [nl] with the WKB approximation, 
give a new exact analytical result, and study the dependence of the spectrum obtained using the 
approximate method of ^ on the choice of the matching time between different time intervals. 
We also comment on recent work subsequent to Q on the trans-Planckian problem for inflationary 
cosmology. 
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I. INTRODUCTION 



The trans-Planckian problem of inflation is the follow- 
ing: in many models of inflation, the phase of accelerated 
expansion lasts so many e-foldings that the comoving 
lengths corresponding to present day cosmological scales 
were much smaller than the Planck length at the begin- 
ning of inflation. Hence, one may wonder whether the 
"standard" predictions of inflation, in particular the fact 
that the power spectrum of cosmological perturbations 
is close to scale-invariant, will be changed if the laws of 
physics beyond the Planck scale are different from the 
ones which rule the low energy phenomena. The usual 
calculations of the spectrum of cosmological perturba- 
tions (see e.g. Q for a comprehensive review) is based on 
the use of classical general relativity coupled to a weakly 
interacting scalar field, and on linearizing the resulting 
equations of motion about a classical background cosmol- 
ogy. The validity of this approach in the trans-Planckian 
regime is highly doubtful. 

In the context of inflation, this question was first ad- 
dressed in Ref. To calculate the power spectrum of 
cosmological perturbations, the main equation that needs 
to be solved is the equation of a parametric oscillator 
with a time-dependent frequency which is a function of 
the scale factor a{rj) (and its derivatives) and of the dis- 
persion relation uj^^^^^{k), where k indicates the physical 
wave number related to the comoving wave number n by 
k = n/a. In the usual discussions, the dispersion relation 
is taken to be linear as appropriate for a free field theory. 
The method used in Ref. was to replace the linear 
relation w . = A; by a non-standard one that mimics 
possible modifications of the physics in the ultraviolet 



regime. Based on similar approaches [^, Q used to study 
the possible dependence of Hawking radiation on trans- 
Planckian physics, two classes of dispersion relations were 
considered [ries = a{'q)uj^^^^{n/ a)\. 
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where A denotes the physical wavelength of a given mode, 
and £^ is a characteristic length expected to be deter- 
mined by the Planck scale. The first one is the Un- 
ruh dispersion relation whereas the second one is the 
Corley/Jacobson relation [Q (6i is an arbitrary number 
which can be positive or negative). In fact, in Ref. a 
generalization of the Corley/Jacobson dispersion relation 
was considered. 

In ||], the problem was investigated for the class of 
scale factors corresponding to power-law inflation, i.e. 
a{ri) — ^ol'yl^^'^i P 1^ ~2 where £q has the dimension 
of a length and is equal to the Hubble radius during in- 
flation if P — —2 (de Sitter inflation). The example of the 
Unruh dispersion relation was treated only for (3 = —2 
whereas the second dispersion relation was studied for 
any value of /3 < — 2. It was found that no modifications 
in the spectrum of fluctuations arise in the first case, 
whereas some differences can show up in the second case 
ff bi < 0\ 

The aim of this letter is to return to the example of 
the Corley/Jacobson dispersion relations, making use of 
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This result was obtained assuming that the initial state is the 
"minimizing energy state". In Ref. another state was also 
considered, but only to demonstrate that the final spectrum de- 
pends on the choice of the initial state. As stressed in Ref. M, the 
"minimizing energy state" is the only physically well-motivated 
state. 
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a new exact analytical solution. This one allows us to ad- 
dress some technical points that have been raised recently 
in the literature and to compare the method of Ref. |] 
with the other methods used in the works subsequent to 
Explicitly, in Section II, we compare the Wigner- 
Kramers-Brillouin (WKB) approximation method used 
in Ref. j5| to calculate the power spectrum with the 
method used in Ref. jl). In section III, we consider 
the exact formula mentioned above which allows us to 
make a smooth transition from the region where the sub- 
Planckian effects are important to the region where the 
dispersion relation is standard. We take advantage of the 
fact that this new solution is also valid for 61 < to dis- 
cuss in more detail than in Ref. [|] the physical meaning 
of complex solutions in the sub-Planckian region. In the 
fourth section, we study the consequences for the match- 
ing time used in the approximate analysis of [Q and show 
that the use of an incorrect matching time will lead to 
artificial oscillations in the spectrum^. First, however, 
we comment on some of the recent work on this subject. 

Following Ref. [Q , there has been a significant amount 
of work on the trans-Planckian problem for inflationary 
cosmology. The work has focused on two issues. The first 
is how broad the class of dispersion relations is for which 
there is a change in the spectrum of fluctuations, and 
whether there are other features in the linear perturba- 
tions which can be used to probe trans-Planckian physics. 
The second is whether the back-reaction of the excess 
fluctuations produced by modifled dispersion relations 
are important and can have interesting consequences. 

Let us for the moment focus on the issue of the class of 
dispersion relations for which a modification of the spec- 
trum of fluctuations is found. In the case of Unruh 
dispersion relations was studied in a de Sitter Universe, 
and no change in the spectrum was found, in agreement 
with the corresponding results in No dispersion re- 
lations of Corley/Jacobson type were studied in l^'^. In 

1^ the trans-Planckian problem of inflationary cosmol- 
ogy was addressed assuming that the mode wave function 
is of WKB-type form. Once again, no changes in the 
spectrum were found. It has been pointed out in |^ that 
for dispersion relations which lead to adiabatic evolution 
of the states on sub-Planckian scales - and WKB states 
fall into this category - there are no changes in the spec- 
trum compared to what is obtained with a linear disper- 
sion relation. As shown in m , the evolution in the case of 



first version of : the incorrect choice of the matching time in 
the Corley/Jacobson case with 61 > led to a spectrum which 
was the usual one times a complicated oscillatory function, in- 
stead of to an unmodified spectrum. 
^ On this basis, the claims of Ref. te| that "contradictory results" 
or "opposite results" were found in Refs. Q and ^ are mislead- 
ing. Indeed, the only common case between these two articles is 
the case of the Unruh dispersion relation in de Sitter spaceti: 
for which exactly the same conclusion was obtained in Refs. 
and 1). 



the Corley-Jacobson dispersion relations with a negative 
sign of bi is not adiabatic, and this is the reason that 
the spectrum is modified. Note that the change in the 
spectrum of fiuctuations in the case of Corley/Jacobson 
dispersion relations with 61 < was confirmed in [Q, 
where these dispersion relations were labeled as "excep- 
tional forms of uj(k)" . Other dispersion relations were 
studied in |, 4 

The question of a possible dependence of the pre- 
dictions for linear cosmological fluctuations on trans- 
Planckian physics was recently considered in the context 
of possible approaches to quantum gravity in |l^, 
H, In ||l|, the starting point was mod- 

iflcations of the commutation relations stemming from 
general considerations of short-distance quantum gravi- 
tational effects, and in |^ consequences of short-distance 
non-commutative geometry was investigated. Interest- 
ing deviations of the spectrum of fluctuations from the 
usual results were found. Speciflcally, the analysis of 
showed that non-Gaussian fluctuations are expected, and 
the analysis of revealed changes to the spectral shape. 

To summarize the current state of knowledge on the 
issue of the possible dependence of the spectrum of fluc- 
tuations on trans-Planckian physics, if the dispersion re- 
lation leads to adiabatic evolution of the vacuum state 
on sub-Planckian length scales, then the spectrum is not 
modifled |^, 1). However, this restriction on the class 
of dispersion relations may exclude the cases of actual 
physical interest String theory, M-theory, non- 

commutative geometry and discrete quantum gravity can 
all lead to much more drastic changes in the effective 
dispersion relation and may hence result in changes in 
the spectrum of ffuctuations which can be probed ob- 
servationally with current and future CMB experiments. 
Recently, two concrete examples where a change in the 
spectrum is obtained were studied in Refs. p^ . 

The second issue raised in recent work, in particular in 
iQ, |l9j is whether, if the power spectrum is indeed modi- 
fled by trans-Planckian physics, there is a back-reaction 
problem ^. The back-reaction problem is the following: 
when viewed at late times, the modified dispersion re- 
lations studied in |^ and in subsequent papers lead to 
mode functions which during inflation on length scales 
larger than £^ but smaller than the Hubble radius are 
excited compared to the adiabatic vacuum. These fluc- 
tuations carry energy and momentum, and this energy 
could be so large as to turn off inffation, in a similar 
manner as the back-reaction of cosmological fluctuations 
in models of chaotic inflation pQ] can build up and ter- 



* The goal of ]lO[ was to nrovide anew model for dark energy (see, 
however, ||l3||), and in and | |l2{ novel dispersion relations were 
^d to yield realizations of the varying speed of light scenario 
and to obtain inflation from radiation [hd]. 
Clearly, the back-reaction problem can only be raised as a second 
question once it has been established that there are important 
modifications to the fluctuation spectrum. 
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minate inflation. This is a very interesting issue which 
merits detailed study. 



II. WKB APPROXIMATION VS. 

SOLUTIONS 



PIECEWISE 



The equation of motion that needs to be solved in order 
to compute the power spectrum of cosmological pertur- 
bations is 



(3) 



where the expression for lo is given below [see Eq. (^]. 
Usually, Eq. (||) possesses two regimes depending on 
whether the wavelength A of the mode is smaller or larger 
than the Hubble radius In the present context, we 
deal with three different regimes. The first one is when 
A < £c < 4- This defines Region I of Ref. §. The sec- 
ond regime is such that <^ \ <^ and corresponds 
to Region II of Ref. §. Finally, the third regime is when 
< 4 < A: this is Region III of Ref. §. We re- 
cover the usual result for the spectrum of fluctuations if 
in Region II /i is given by a single branch with a coef- 
ficient proportional to \fn multiplying the plane wave 
solution. 

The method used in Ref. Q to determine the solution 



of Eq. (||) is to find solutions in the three regions and to 
match them at the boundaries between the regions. The 
method used in Ref. ||^ is different and consists in utiliz- 
ing the WKB approximation. In that case, the solution 
of Eq. (||) is given by 



Mwkb 



exp 



±i / w(r)dT 



(4) 



where r\\ is some initial time. 

In cases where the adiabatic approximation is justified, 
the two methods are equivalent. To verify this, consider 
first Unruh's dispersion relation. Inserting the disper- 
sion relation into Eq. (^) (using the minus sign) yields in 
Region I 



-27ri/£-|-l/2 



(5) 



This is nothing but Eq. (66) of Ref. || [i.e. Eq. (54) in 
the limit when t = 1,^1 vo small], properly normalized 
according to Eqs. (56)- (57) of that reference. Similarly, 
in the case of the Corley/Jacobson dispersion relation 
(with h\ > 0, in which case the adiabatic condition is 
satisfied [O), Eq. (0) becomes 



Mwkb 



(6) 



Again, this is the same result as obtained with the match- 
ing technique in Ref. [0|. To see this, take Eq. (118) of 
that reference, which provides the solution in Region I, 
consider the limit when z is large and use Eqs. (126) and 
(127) with the lower sign. In Region II, the solution is 
simply given by plane waves /i ~ B\e™'^ -I- B2e^™^ in 
agreement with Eq. (14) of Ref. M. 



The advantage of the WKB method over the method 
used in Ref. [|l[ is that one does not have to perform the 
matching between the different regions. The disadvan- 
tage is that the WKB method applies only to examples 
in which the adiabaticity condition is satisfled. In the 
next section, we consider an exact model which allows us 
to avoid the matching between Region I and II. In the 
third section, we study the consequences for the match- 
ing time and show that an incorrect matching time could 
lead to artificial oscillations in the spectrum, as pointed 
out in Ref. g. 



III. THE EXACT MODEL 

The effective time dependent frequency is given by the 
equation 



a 
a 



(7) 



We consider the Corley-Jacobson dispersion relation 
given in Eq. (|^) where, as already mentioned above, bi is 
an arbitrary number which can be positive or negative. 
We restrict our consideration to the prototypical model 
of inflation, i.e. de Sitter inflation. Then the scale factor 
can be written as 0(77) — £Q/\ri\. If we consider the evolu- 
tion of the mode well inside horizon, the term a" /a can 
be neglected in Eq. (^). Then, this equation takes on the 
form 



2 \ 2 



(8) 



where e = i^/io- Typically, e is a very small number of 
the order 10~^. We now need to distinguish between the 
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cases of positive or negative 61. Let us first concentrate 
on the 61 > case. We make the following change of 
variables — ry = Cx (a; > 0), where C is a constant given 
by C = [{bi)~'^ ^ ^ / n] ■\/7r/e. Then the equation of motion 
takes the form 

Clearly, this equation possesses two regimes. The first 
one, corresponding to Region I in Ref. [Q, is when the 
quartic term of the dispersion relation dominates. Ex- 
pressed in terms of the new variable x, this corresponds 



to x^/4 > |a|. On the other hand, Region II of Ref. @] 
corresponds to a region where the dispersion relation has 
become standard and a;-^/4 <C \a\. The general solution 
of Eq. (^) (valid in both regions) is given in terms of 
parabohc cyhnder functions E(a,x) (see Ref. |^) and 
reads 

AiE{a,x) + A2E*{a,x), (10) 

where the two constants Ai and A2 are fixed by the initial 
conditions when 77 — > —00 (i.e. x — > 00). In this limit, 
the asymptotic behavior is given by 



lim /X 

rf — > — 00 



-1/8 



+A2 cxp 



-1/4 



-1/2 



H-'/Ha, 



_ exp 



•1,1/2"- e, r 



I I — 

4 2 



where (f)2 = argr(l/2 -f ia). Thus we reproduce the correct WKB behavior if we take 



, ,,-1/8 /^X 1/4 



exp 



4 2 



^2=0. 



r 



(11) 



(12) 



Let us now study how this solution behaves in Re- 
gion II. Usually, the solution in this region is given by 
l/\/2ne^™l''l. It is therefore sufficient to have only one 
branch proportional to 1/Vn to recover the standard 
scale-invariant spectrum. Using the asymptotic behav- 
ior of the parabolic cylinder functions when x^ /A<^ |a|, 
one finds 



-1/4 



1/8 



exp 



in\ri\ + i— 



(13) 



i.e. we precisely recover the conditions necessary to ob- 
tain a Harrison- Zeldovich spectrum, due to the fact that 
Ax cx l/^/n, see Eq. (^. 

Let us now turn to the case 61 < 0. An immediate con- 
sequence is that the dispersion relation rieff (n, r]) vanishes 
at some point and then becomes complex. Therefore, we 
face the following problems. Firstly, in the region where 
nee{n,r]) is small, the term a" /a is no longer neghgible. 
In principle it should be taken into account in the equa- 
tion of motion but no exact solution can then be found. 
However, this problem is not too serious because the ef- 
fect of this term on the final spectrum is expected to be 
small. Secondly, a much more serious question is the fact 
that there is a region where we have to quantize a field 
in the presence of imaginary frequency modes. Although 
imaginary frequencies are standard in classical physics 
and in quantum mechanics, they seem to be problematic 
in the context of quantum field theory |23] although there 



exist concrete physical situations where they are impor- 
tant . Thirdly, we have to fix the initial conditions in 
the complex region. A possible choice for the initial con- 
ditions, which is consistent with the WKB result, is to 
keep only the decreasing exponential in the region where 
the effective frequency becomes complex as proposed in 
Ref. [|] . Then, an exact solution in terms of parabolic 
cyhnder functions U{a,x) and V{a,x) (see Ref. [^) can 
be found. This gives a Harrison- Zeldovich spectrum cor- 
rected by oscillations and by an exponential term of the 
form e"^" . This result is in agreement with what was 
obtained in Ref. However, this result rests clearly 
on "non-standard physics" and for this reason is not so 
attractive. Therefore, it is important to recall that there 
now exist two cases where the final spectrum is modified 
and everything is well-defined lisl . 



IV. CONSEQUENCES FOR THE MATCHING 
PROCEDURE 



We now investigate what we can learn from the previ- 
ous exact solution with respect to the matching between 
Regions I and II. We concentrate on the Corley/Jacobson 
case with 61 > 0. A priori, two natural choices for the 
matching time can be envisaged. The first choice is to 
match the solutions when X = £„■ This amounts to 
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choosing 

Another possibihty is to choose the matching time such 
that the usual and the extra contribution in the disper- 
sion relation are equal 




(15) 



This is in fact equivalent to matching the frequencies 
nosin,-!]). Therefore, 771 and ry^ are not equal unless &i = 
1. Let us now perform the matching between Regions I 
and II. Then, the coefficient Bi is given by: 



y -1/2 

^<^^ -1/2,-1/4, 1-1/2 



exp 



1/2/ ne 



x|^n-l(l + /3)hJrl±^6}/^(^g)n|r;j| 



For the curly bracket of the above expression one finds 



{•■■}i = -^(l+/3)l^jr'+*n(l±|||) 



(16) 



(17) 



In the same manner, one can determine B2 to be 



X < m 



m+i(l+/3)|,yjriT*6y'(g)nhl}. 



(18) 



Again, we can write the curly bracket in the previous 
expression as 



{•••}2 



1 



m 1 =F 



(19) 



The situation is now clear. If the joining is performed at 
rij =771, the terms proportional to in in the curly brack- 
ets have no reason to cancel out and they are in fact of 
order one. The term proportional to |7?j|~^ is very small 
and can be neglected. Therefore we reach the conclu- 
sion that El ~ B2 and we have oscillations. However, 
if we perform the matching at r/j = ri[ the situation is 
drastically different. This time one of the curly brack- 
ets vanishes and one of the i?i's becomes of the order 
rj-j^^ <C 1, whereas the other one is of order 1. In other 
words, only one branch survives, we have no oscillations 
and the spectrum of fiuctuations is unchanged. 

A comparison with the exact solution shows that the 
correct matching time is rj'i- Therefore, there are no oscil- 
lations in the spectrum in the case 61 > 0, as pointed out 
in Ref . |^ , and the spectrum is unchanged. On the other 
hand, for bi < the spectrum is modified, in agreement 



with the analysis of Ref. |^ . 



V. CONCLUSION 

In this short letter, we have investigated the following 
technical points: (i) it has been shown that the method 
used in Ref. is equivalent to the WKB approach, (ii) 
a new solution valid in the case of the Corley/Jacobson 
dispersion relation has been presented, (iii) in the case 
where piecewise solutions are used, the matching condi- 
tions have been studied. It has been demonstrated that 
the frequencies rather than the wavelengths should be 
matched (as could have been guessed from the equation of 
motion) and that if the latter requirement is utilized then 
artificial oscillations show up in the spectrum. These re- 
sults complete the study of Ref. |jl| . 

On more general grounds, the conclusion that follows 
from the previous considerations is that there is a sen- 
sitive dependence of the spectrum of cosmological fluc- 
tuations on the assumptions made at the level of sub- 
Planckian physics. A separate issue is whether these 



6 



modifications are reasonable from a physical point of 
view. This question cannot be answered in the absence 
of a realistic theory of physics beyond the Planck scale. 
To our knowledge, it is still an open problem to derive a 
dispersion relation from, for example, string theory (see, 
however, Refs. ^ for some recent progress). 

Acknowledgements 

We are grateful to Steve Corley for stimulating dis- 



cussions and useful comments. We acknowledge support 
from the BROWN-CNRS University Accord which made 
possible the visit of J. M. to Brown during which some 
of the work on this project was done, and we are grate- 
ful to Herb Fried for his efforts to secure this Accord. 
J. M. thanks the High Energy Group of Brown Univer- 
sity for warm hospitality. The research was supported in 
part by the U.S. Department of Energy under Contract 
DE-FG02-91ER40688, TASK A. 



[1] 



[2] 

[3] 
[4] 



J. Martin and R. H. Brandenbe iRer, Phys. Rev. D 63, 
123501 (2001), |iep-th/0005209|; R. H. Brandenberger 



and J. Martin, M od. Phys. Lett. A 16, 999 (2001), 
astro- ph/0005434 



V. F. Mukhanov, H. A. Foldman and R. H. Branden- 
berger, Phys. Rept. 215, 203 (1992). 
W. Unruh, Phys. Rev. D 51, 2827 (1995). 
S. Corley and T. Jacobson, Phys. Rev. D 54, 1568 (1996), 
hep-th/9601073|. 



J. C. Niemey er and R. Parenta ni, Phys. Rev. D 64, 
101301 (2001), [Iit"ro-ph/0101451 . 

J. C. Niemeyer, Phys. Rev. D 63, 123502 (2001), 
astro- ph/0005534 



[8] 
[9] 

[10] 

[11] 
[12] 

[13] 



(2001) [JETP Lett. 73, 371 (2001)], astro-ph/0104043 . 
J. Martin and R. H. Brandenberger, astro-ph/0012031 . 
J. Kowalski-Glikni an, Phys. Lett. B 499, 1 (2001), 
astro- ph/000625c| . 



L. Mersini, M. Bas tero-Gil and P. K anti, Phys. Rev. D 
64, 043508 (2001), ^i^-ph/010121o|. 



S. Alexander and J. Magueijo, |hep-th/010409; . 
S. Alexander. R. Brandenberger and J. Magueijo, 



hep-th/0108190 . 

M. Lemoine, M. Lubo, J. Martin and J. Uzan, Phys. Rev. 



[14] 

[15] 

[16] 

[17] 

[18] 

[19] 
[20] 

[21] 



[22] 
[23] 

[24] 



D 64, 043508 (2001), ^iep-th/0109128| . 
A. Kompf, Phy s. Rev. D 63, 083514 (2001), 



astro-ph/0009209 



C. Chu, B. R. G reene and G. Shi u, Mod. Phys. Lett. A 



16, 2231 (2001), hep-th/0011241 



A. Ke mpf and J. C. Nieme yer, Phys. Rev. D 64 103501 
(2001), ^stro-ph/0103224 

R. Easther, B. R. Greene, W. H. Kinney and G . Shiu, 
Phys. Rev. D 64 103502 (2001), ^ep-th/0104104 
R. H. Branden berger, S. E. Joras and J. Martin, 



hep-th/0 112122 



astro-ph/0012431 



T. Tanaka, 

L. R. W. Abramo, R. H. Brandenberger and 
V. Mukhanov, Phys. Rev. D 56, 3248 (1997); 
V. Mukhanov, L. R. W. Abramo and R. H. Branden- 
berger, Phys. Rev. Lett. 78, 1624 (1997). 
I. S. Gradshteyn, I. M. Ryzhik, "Table of Integrals, 
Series, and Products", Sixth edition (2000), Academic 
Press. 

G. Gang, liep-ph/9603166 . 

L. I. Schiff, H. Snyder and J. Weinberg, Phys. Rev. 57 
315 (1940) 

G. Amelino-Camelia, gr-qc/0012051 



